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Abstract 

We describe the sequences {x n } n given by the non-autonomous second order Lyness 
difference equations x n+ 2 = (a„ + x n+ i)/x n , where {a„}„ is either a 2-periodic or a 3- 
periodic sequence of positive values and the initial conditions x%, xi are as well positive. 
We also show an interesting phenomenon of the discrete dynamical systems associated 
to some of these difference equations: the existence of one oscillation of their associated 
rotation number functions. This behavior does not appear for the autonomous Lyness 
difference equations. 
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1 Introduction and main result 

This paper fully describes the sequences given by the non-autonomous second order Lyness 
difference equations 

(In , - 

X n +2 = , (1) 

where {a n } n is a fc-periodic sequence taking positive values, k = 2,3, and the initial con- 
ditions x\,X2 are as well positive. This question is proposed in [H Sec. 5.43]. Recall 
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that non-autonomous recurrences appear for instance as population models with a variable 
structure affected by some seasonality [IT] , where k is the number of seasons. Some 
dynamical issues of similar type of equations have been studied in several recent papers 

m Big nana us ma. 

Recall that when k = 1, that is a n = a > 0, for all n £ N, then (JTJ) is the famous Lyness 
recurrence which is well understood, see for instance [21 [18]. The cases k = 2, 3 have been 
already studied and some partial results are established. For both cases it is known that 
the solutions are persistent near a given /c-periodic solution, which is stable. This is proved 
by using some known invariants, see [14 \ [TT] . Recall that in our context it is said that a 
solution {x n } n is persistent if there exist two real positive constants c and C, which depend 
on the initial conditions, such that for &U n > 1,0 < c < x n < C < oo. We prove: 

Theorem 1. Let {x n } n be any sequence defined by and k 6 {2, 3}. Then it is persistent. 
Furthermore, either 

(a) the sequence {x n } n is periodic, with period a multiple of k; or 

(b) the sequence {x n } n densely fills one or two (resp. one, two or three) disjoint intervals 
o/M + when {a n } n is 2-periodic (resp. 3-periodic) . Moreover it is possible by algebraic 
tools to distinguish which is the situation. 

Our approach to describe the sequences {x n } n is based on the study of the natural 
dynamical system associated to ([I]) and on the results of [6]. The main tool that allows 
to distinguish the number of intervals for the adherence of the sequences {x n } n is the 
computation of several resultants, see Section [H 

It is worth to comment that Theorem [1] is an extension of what happens in the classical 
case k = 1. There, the same result holds but in statement (b) only appears one interval. Our 
second main result will prove that there are other more significative differences between the 
case k = 1 and the cases k = 2,3. These differences are related with the lack of monotonicity 
of certain rotation number functions associated to the dynamical systems given by the 
Lyness recurrences, see Theorem [3l The behaviors of these rotation number functions are 
important for the understanding of the recurrences, because they give the possible periods 
for them, see [21 El EH] . 

On the other hand in Q2] it is proved that, at least for some values of {a n } n , the 
behaviour of {x n } n for the case k = 5 is totally different. In particular unbounded positive 
solutions appear. In the forthcoming paper [7] we explore in more detail the differences 
between the cases k = 1,2,3 and k > 4. 

This paper is organized as follows: Section [2] presents the difference equations that 
we are studying as discrete dynamical systems and we state our main results on them, see 
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Theorems [2] and [3l Section [3] is devoted to the proof of Theorem [2j By using it, in Section SI 
we prove Theorem Q] and we give some examples of how to apply it to determine the number 
of closed intervals of the adherence of {x n } n . In Section [5] we demonstrate Theorem [3] and 
we also present some examples where we study in more detail the rotation number function 
of the dynamical systems associated to (HJ. 



2 Main results from the dynamical systems point of view 

In this section we reduce the study of the sequence {x n } n to the study of some discrete 
dynamical systems and we state our main results on them. 
First we introduce some notations. When k = 2, set 



a for n = 21 + 1, 
b for n = 2£, 



and when k = 3, set 



(2) 



a for n = 31 + 1, 

b for n = 3£ + 2, (3) 
c for n = 3£, 

where £ G N and a > 0, b > and c > 0. 

We also consider the maps F a (x, y), with a £ {a, b, c}, as 

F a (x,y) = ( y, ° ^ U 



x 

defined on the open invariant set Q + := {(x, y) : x > 0, y > 0} C M 2 . 
Consider for instance k = 2. The sequence given by (H|), 

x 1 ,x 2 ,x 3 , x 4 ,x 5 , x 6 , x 7 , . . . , (4) 

can be seen as 

(xi,x 2 ) (X2,X 3 ) (£3, x 4 ) (x 4 , x 5 ) (x 5 , x 6 ) • • • . 

Hence the behavior of ([U) can be obtained from the study of the dynamical system defined 
in Q + by the map: 

' a + y a + bx + y " 



Fb,a(x,y) ■= F b oF a (x,y) 
Similarly, for k = 3 we can consider the map: 



x xy 



. a + bx + y a + bx + y + cxy 

F c ,b,a(Xi V) ■= F c oF b o F a (x, y) = , — ■ — r 

xy y(a + y) 
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Notice that both maps have an only fixed point in Q + , which depends on a, b (and c), that 
for short we denote by p. 

It is easy to interpret the invariants for (pQ) and k = 2,3, given in |12l I14j . in terms of 
first integrals of the above maps, see also Lemma [6j We have that 

rr . . ax 2 y + bxy 2 + bx 2 + ay 2 + (b 2 + a)x + (b + a 2 )y + ab 

V b) a{x,y) := , 

xy 

is a first integral for F},^ and 

. . cx 2 y + axy 2 + bx 2 + by 2 + (a + frc)x + (c + a6)« + ac 

xy 

is a first integral for F c & ja . The topology of the level sets of these integrals in Q + as well as 
the dynamics of the maps restricted to them is described by the following result, that will 
be proved in Section [3l 

Theorem 2. (i) The level sets of Vb, a (resp. V Cj b,a ) ^ n Q + \ {p} are diffeomorphic to 
circles surrounding p, which is the unique fixed point of F^a (resp. F c> b a ). 

(ii) The action of F\, a (resp. F c ^, a ) on each level set of Vb. a (resp. V c ,b.a) contained in 
Q + \ {p} is conjugated to a rotation of the circle. 

Once a result like the above one is established the study of the possible periods of 
the sequences {x n } n given by ([T]) is quite standard. It suffices, first to get the rotation 
interval, which is the open interval formed by all the rotation numbers given by the above 
theorem, varying the level sets of the first integrals. Afterwards, it suffices to find which are 
the denominators of all the irreducible rational numbers that belong to the corresponding 
interval, see [31 El [18] . 

The study of the rotation number of these kind of rational maps is not an easy task, 
see again [21 [31 118] . In particular, in [2] was proved that the rotation number function 
parameterized by the energy levels of the Lyness map F a ,a ^ 1, is always monotonous, 
solving a conjecture of Zeeman given in [18], see also [15]. As far as we know, in this paper 
we give the first simple example for which this rotation number function is neither constant 
nor monotonous. We prove: 

Theorem 3. There are positive values of a and b, such that the rotation number function 
Pb,a(h) of Fb t a associated to the closed ovals of {Vb A = h} C Q + has a local maximum. 

Hence, apart from the known behaviors for the autonomous Lyness maps, that is global 
periodicity or monotonicity of the rotation number function, which trivially holds for Ff,^, 
taking for instance o = 6= lora = 6/l, respectively, there appear more complicated 
behaviors for the rotation number function. 
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Our proof of this result relies on the study of lower and upper bounds for the rotation 
number of F^ a on a given oval of a level set of V^ a given for some (a, b) G (Q + ) 2 and 
{Vb,a{x,y) = Vb,a(%o, 2/o)}> for (xo,yo) S (Q + ) 2 . This can be done because the map on this 
oval is conjugated to a rotation and it is possible to use an algebraic manipulator to follow 
and to order a finite number iterates on it, which are also given by points with rational 
coordinates. So, only exact arithmetic is used. A similar study could be done for F C ^ A . 

3 Proof of Theorem [2] 

Proof of (i) of Theorem^ The orbits of Fb ia and F C ^ A lie on the level sets Vb t a = h and 
V Cj b,a = h respectively. These level sets can be seen as the algebraic curves given by 

C*2 := { c 2{x, y) = ax 2 y + bxy 2 + bx 2 — hxy + ay 2 + (b 2 + a)x + (b + a 2 )y + ab = 0} 

and 

C?j := {c3(x, y) = cx 2 y + axy 2 + bx 2 — hxy + by 2 + (a + bc)x + (c + ab)y + ac = 0}, 
respectively. 

Taking homogeneous coordinates on the projective plane PM. 2 both curves C2 and C3 
have the form 

C := {Sx 2 y + Txy 2 + Ux 2 z + Vxyz + Wy 2 z + Lxz 2 + Myz 2 + Nz 3 = 0}. 

In order to find the branches of them tending to infinity, we examine the directions of 
approach to infinity {z = 0) in the local charts determined by x = 1 and y = 1 respectively. 
In the local chart given by x = 1, the curve C writes as 

Sy + Ty 2 + Uz + Vyz + Wy 2 z + Lz 2 + Myz 2 + Nz 3 = 

and it meets the straight line at infinity z = when y(S + Ty) = 0. Since for both curves 
Ci and C3 the coefficients S and T are positive, the only intersection point that could give 
points in Q + is (y,z) = (0,0). The algebraic curve C arrives to (y,z) = (0,0) tangentially 
to the line Sy + Uz = 0. Since for both curves, C2 and C3, the coefficients S and U are 
also positive, we have that the branches of the level sets tending to infinity are not included 
in Q+. 

An analogous study can be made in the chart given by y = 1, obtaining the same 
conclusions. 

Moreover, it can be easily checked that in the affine plane both curves C2 and C3 do 
not intersect the part of the axes x = and y = which is in the boundary of Q + . 
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In summary, there are no branches of the curves C2 and C3 tending to infinity or crossing 
the axes x = and y = in Q + , and therefore the connected components of Cj n Q + for 
i = 2, 3 are bounded. Notice that this result in particular already implies the persistence of 
the sequences given by ([TJ. 

Consider k = 2. We claim the following facts: 

(a) In Q + , the set of fixed points of Fb <a and the set of singular points of C% coincide and 
they only contain the point p = 

(b) The function Vb t a(x,y), has a local minimum at p. 

We remark that item (b) is already known. We present a new simple proof for the sake of 
completeness. 

From the above claims and the fact that the connected components of the level sets of 
Vb A in Q + are bounded it follows that the level sets of Vb >a in Q + \ {p} are diffeomorphic 
to circles. 

Let us prove the above claims. The fixed points of Ff, <a are given by 
x = 2±K 

X ' 




£t£±at, \x{y 2 -b) = a + y 



xy 

and so x 2 = x(y 2 — b). Hence in Q + , we have that x = y 2 — b and the above system is 
equivalent to 

x = y 2 — b, I x = y 2 — b, 

xy 2 - bx - y - a = 0, I P(y) := y 4 - Iby 2 - y + b 2 - a = 0. 

It is not difficult to check that the last system of equations is precisely the one that gives 
the critical points of the curves V a b = h. Moreover, from the first equation it is necessary 
that x = y 2 — b > and hence y > y/b. Since P(y) has only one real root in (Vb, 00) the 
uniqueness of the critical point holds. 

Let us prove that this critical point corresponds with a local minimum of Vf, a . We will 
check the usual sufficient conditions given by the Hessian of at p. 

Firstly, 

^ tt (2 t \ (y + a) {ay + b) / - 
— 2 V b , a (y -b,y) = 2 ^_ hfy >0for y >A 

Secondly, the determinant of the Hessian matrix at the points (y 2 — b, y) is 

(b - y l Yy^ 



6 



where 



f{y) := {by 2 + a - b 2 ){-by 6 + 3 (a + b 2 ) y 4 + 4 (a 2 + b) y 3 + 36 (2a - 6 2 ) y 2 + b 2 (b 2 - a). 



r(y) = (4a 2 6 2 + 4a 3 + 46 3 + 6a&) y 3 + (l8a 2 6 + 86 3 a + 36 2 ) y 2 

+ (-46 3 a 2 + 4a 3 6 - 46 4 + 12ao 2 + 3a 2 ) y - 86 4 a + 5a 2 b 2 + 3a 3 . 

Observe that if y is the positive root of P(y), then sign(/t(y)) = sign(r(y)). Taking into 



r(y) = y 2 (4y 3 - 46y - l) (by + 1 - y 3 ) 2 (6 - ff . 
So, sign (4y 3 - 4by - l) = sign (P'(y)) . Since P(v / 6) = - a < 0, lim P(y) = +oo 



and, on this interval, there is only one critical point of P(y), which is simple, we get that 
P'(y) > and so h(y) > 0. Hence p is a local minimum of Vft ja (x, y), as we wanted to prove. 

The same kind of arguments work to end the proof for the case k = 3, but the compu- 
tations are extremely more tedious. We only make some comments. 

The fixed points of F c ^ a in Q + are given by: 



It can be proved again that they coincide with the singular points of V c \,,a i n Q + ■ This fact 
follows from the computation of several suitable resultants between dV c ^ a /dx, dV c j> %a /dy 
and V c ,b, a - 

The uniqueness of the fixed point p in Q + can be shown as follows: since Q(x) > 
implies that x > 1, we only need to search solutions of P^{x) = in (1, +oo). With the new 
variable z = x — 1, 



P(z) := P(z + 1) := z 5 + (c + 5)z 4 + (8 + 4c)z 3 + (4-ab + 4c)z 2 - (a + b)z - (a + b) 2 = 0, 



Since P(0) < 0; lim P(z) = +oo; and the Descarte's rule, we know that there is only one 



positive solution, as we wanted to see. 

Finally it can be proved that p is a non-degenerated local minimum of V c f, a . These 
computations are complicated, and they have been performed in a very smart way in |14j . 



A tedious computation shows that f(y) = q{y)P(y) + r(y), with 



account that P{y) = implies that a = y A — 2by 2 — y + b 2 



we have that 




■ 4 - 2x 3 - (2a + bc)x 2 + (1 - b 2 - c 2 )x + a - be = 0, 



so we skip them and we refer the reader to this last reference. 



I 
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3.1 Proof of (ii) of Theorem [2] 

In [6] it is proved a result that characterizes the dynamics of integrable diffeomorphisms 
having a Lie Symmetry, that is a vector field X such that X(F(p)) = {DF(p)) X(p). Next 
theorem states it, particularized to the case we are interested. 

Theorem 4 ([B]). Let WcM 2 be an open set and let 3> :U —> U be a diffeomorphism such 



(a) It has a smooth regular first integral V : U — * M, having its level sets I\ := {z = 
(x, y) £U : V{z) = h} as simple closed curves. 

(b) There exists a smooth function /i :U —* M + such that for any z £W, 



Then the map restricted to each I\ is conjugated to a rotation with rotation number 
r(h) /T(h), where T(h) is the period ofT^ as a periodic orbit of the planar differential 
equation 



and t(/i) is the time needed by the flow of this equation for going from any w G to 



Next lemma is one of the key points for finding a Lie symmetry for families of periodic 
maps, like the 2 and 3-periodic Lyness maps. 

Lemma 5. Let {G a } a ^A be a family of diffeomorphisms of U C M?. Suppose that there 
exists a smooth map ^ : U — > M such that for any a £ A and any z £ U, the equation 
fi(G a (z)) = det(DG a (z)) /j,(z) is satisfied. Then, for every choice a%, . . . , a& G A, we have 



where = G ak o ■ ■ ■ o G a2 o G ai . 

Proof. It is only necessary to prove the result for k = 2 because the general case follows 
easily by induction. Consider ai,a2 G A then 



/i(G a2 , 0l (z)) =n(G a2 °G ai (z)) = det(DG a2 (G ai (z)))fi(G ai (z)) = 

= det(DG a2 (G ai (z))) det(DG ai (z)) n(z) = det(D(G a2 o G ai (z))) fj,(z) = 
= det(DG a2 , ai {z))n(z), 



that: 



fi($(z)) = det(D$(z))n(z). 




&(w) G r h . 



K G [k](z)) = det(DG[ k ](z)) fi(z) 



and the lemma follows. 



I 
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Proof of (ii) of Theorem O From part (i) of the theorem we know that the level sets 
of V& )Q and V c>b>a m Q + \ {p} ar e diffeomorphic to circles. Moreover these functions are 
first integrals of F bi0j and F Cjbja , respectively. Notice also that for any a, the Lyness map 
F a (x,y) = (y, ^±J/) satisfies 

K F a(x,y)) = dek(DF a (x,y))fi(x,y), 
with n(x,y) = xy. Hence, by Lemma [5j 
fJ,(F b>a (x, y)) = &et{DF b)CL (x, y))n(x, y) and n(F c ^ a (x, y)) = &et{DF c ^ a (x, y))fJt(x, y). 

Thus, from Theorem 01 the result follows. g 

It is worth to comment that once part (i) of the theorem is proved it is also possible to 
prove that the dynamics of F b>a (resp. F C)b)a ) restricted to the level sets of V bjCL (resp. V CjbA ) 
is conjugated to a rotation by using that they are given by cubic curves and that the map is 
birational, see [T3]. We prefer our approach because it provides a dynamical interpretation 
of the rotation number together with its analytic characterization. 

4 Proof of Theorem [L] 

In order to prove Theorem [1] we need a preliminary result. Consider the maps Fb^ a and 
F a fi, jointly with their corresponding first integrals Vb t a and V ab - in a similar way consider 
F c ,b,a , F a ^ h and F b)a>c with V Cjbja , V a ^ b and V b ^ c . Some simple computations prove the 
following elementary but useful lemma. Notice that it can be interpreted as the relation 
between the first integrals and the non-autonomous invariants. 

Lemma 6. With the above notations: 

(i) Vb,a( X >y) = V a,b{F a {x,y)). 

(H) Vc !bt a(x,y) =V a)C:b (F a (x,y)) =V b , a , c (F b (F a (x,y))). 

Proof of TheoremUl We split the proof in two steps. For k = 2,3 we first prove that there 
are only two types of behaviors for {x n } n , either this set of points is formed by kp points 
for some positive integer p, or it has infinitely many points whose adherence is given by at 
most k intervals. Secondly, in this later case, we provide an algebraic way for studying the 
actual number of intervals. 

First step: We start with the case k = 2. With the notation introduced in ([2]), it holds 
that 

F bt a{x2n-l,X2n) = (^2n+l , X2 n +2) , F ajb (x2 n , »2n+l) = { x 2n+2, X2n+3), 
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where (xi,^) G Q + and n > 1. So the odd terms of the sequence {x n } n are contained in 
the projection on the x-axis of the oval of {Vft ja (x, y) = Vj, j0 (xi, X2) = h} and the even terms 
in the corresponding projection of {V a) b(x, y) = V a ^{F a {x\, X2)) = /i}, where notice that we 
have used Lemma [6l 

Recall that the ovals of Vb, a are invariant by F ba and the ovals of V a f> are invariant by 
-F a ,b- Notice also that the trivial equality F a o a = F a j, o F a implies that the action of F;, ja 
on {Vft ia (x,y) = /i} is conjugated to the action of F a> b on {^^(x, = /i} via F a . 

From Theorem[2]we know that F b a on the corresponding oval is conjugated to a rotation 
of the circle. Hence, if the corresponding rotation number is rational, then the orbit starting 
at (zi, X2) is periodic, of period say q, then the sequence {x n } n is 2c/-periodic. On the other 
hand if the rotation number is irrational, then the orbit of (xi, X2) generated by F b a fulfills 
densely the oval of {Vb A (x,y) = h} in Q + and hence the subsequence of odd terms also 
fulfills densely the projection of {Vb t a{x, y) = h} in the x-axis. Clearly, the sequence of even 
terms do the same with the projection of the oval of {V a b(x,y) = h}. 

Similarly when k = 3 the equalities 



where n > 1, allow to conclude that each term x m , of the sequence {x n } n where we 
use the notation ([3]), is contained in one of the projections on the x-axis of the ovals 
{V c ,b, a (x,y) = V C)bta (xi,X2) =■■ h} and {V a ^ b {x,y) = h} and {V b , a , c (x,y) = h}, according 
with the remainder of m after dividing it by 3. The rest of the proof in this case follows as 
in the case k = 2. So the first step is done. 

Second step: From the above results it is clear that the problem of knowing the number 
of connected components of the adherence of {x n } n is equivalent to the control of the pro- 
jections of several invariant ovals on the x-axis. The strategy for k = 3, and analogously for 
k = 2, is the following. Consider the ovals contained in the level sets given by {V Ct b, a (x, y) = 
Mi {Ya,c,b{ x iV) = h} an d {yb,a,c{ x iV) = M an d denote by I = I(a,b,c,h), J = J(a,b,c,h) 
and K = K(a, b, c, h) the corresponding closed intervals of (0, 00) given by their projections 
on the x-axis. 

We want to detect the values of h for which two of the intervals, among /, J and K, have 
exactly one common point. First we seek their boundaries. Since the level sets are given by 
cubic curves, that are quadratic with respect the y-variable, these points will correspond 
with values of x for which the discriminant of the quadratic equation with respect to y is 



A c,b,a n—2i XSn—1 ) 
FajC,b\ x 3n— 1 ? ^3ra) 
± b,a,c n > 2-3n+l J 



{Xs n+ i, X3 n+ 2) 
(x3n+2,X3 n+ 3) 
(^371+3, X3 n+ 4) 
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zero. So, we compute 



Ri(x, h, a, b, c) := dis (xyV c ^ a (x, y) - hxy, y) = 0, 
R 2 (x, h, a, b, c) := dis (xyV a ^ b (x, y) - hxy, y) = 0, 
R 3 (x, h, a, b, c) := dis {xyV b ^ c {x, y) - hxy, y) = 0. 

Now we have to search for relations among a, b, c and h for which two of these three functions 
have some common solution, x. These relations can be obtained by computing some suitable 
resultants. 

Taking the resultants of Ri and R 2 ; R 2 and R%; and R± and R% with respect x we 
obtain three polynomial equations R^(h, a, b, c) = 0, Rs(h, a, b, c) = and Re(h, a, b, c) = 0. 
In short, once a, b and c are fixed we have obtained three polynomials in h such that a 
subset of their zeroes give the bifurcation values which separate the number of intervals 
of the adherence of {x n } n . See the results of Proposition [7| and Example [5] for concrete 
applications of the method. 

Before ending the proof we want to comment that for most values of a, b and c, varying 
h there appear the three possibilities, namely 1, 2 or 3 different intervals. The last case 
appears for values of h near h c := V Ct b, a (p), because the first coordinates of the three points 
p, -F a (p) and F{,(F a (p)) almost never coincide. The other situations can be obtained by 
increasing h. g 

Proposition 7. Consider the recurrence (pQ) with k = 2 and {a n } n as in (j2]) taking the 
values a = 3 and b = 1/2. Define h c = (12z 3 - 33z + 7)/(2(z 2 - 3)) ~ 17.0394, where 
z ~ 2.1513 is the biggest positive real root of 2z 4 — 12z 2 — 2z + 17, and h* ~ 17.1198, is 
the smallest positive root of 

p 4 (h) := 112900/i 4 - 2548088/i 3 - 48390204/i 2 + 564028596/i + 7613699255. 

Then, 

(i) The initial condition {x\,x 2 ) = (z,z 2 — 3) gives a two periodic recurrence {x n } n . 
Moreover ^1/2,3(2, 3 — z 2 ) = h c . 

(ii) Let (x±,x 2 ) be any positive initial conditions, different from (z,z 2 — 3), and set h = 
Vi/2,3 ( x i j x 2)- Let p(h) denote the rotation number of -Pi/2,3 restricted to the oval of 
{Vi/ 2 , 3 {%, y) = h}. Then 

(I) If p(h) =p/q £ Q, with gcd(p,q) = 1, then the sequence {x n } n is 2q-periodic. 

(II) If p(h) Q and h E (h c ,h*) then the adherence of the sequence {x n } n is formed 
by two disjoint closed intervals. 
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(Ill) If p(h) Q and h £ [h*,oo) then the adherence of the sequence {x n } n is one 
closed interval. 

We want to remark that, from a computational point of view, the case (I) almost never 
is detected. Indeed, taking a and 6 rational numbers and starting with rational initial 
conditions, by using Mazur's theorem it can be seen that almost never the rotation number 
will be rational, see the proof of [SJ Prop. 1]. Therefore, in numeric simulations only 
situations (II) and (III) appear, and the value h = h* gives the boundary between them. 
In general, for k = 2, the value h* is always the root of a polynomial of degree four, which 
is constructed from the values of a and b. 

Proof of Proposition^ Clearly (z, 3—z 2 ) is the fixed point of F^ a in Q + . Some computations 
give the compact expression of h c := V a j,{z, z 2 — 3). To obtain the values h* we proceed as 
in the proof of Theorem [TJ In general, 

Ri(x, h, a, b) := dis (xyV b , a (x, y) - hxy, y) 

= (ax 2 -hx + a 2 + b) 2 - 4(6x + a)(bx 2 + b 2 x + ax + ab), 
R 2 (x, h, a, b) := dis (xyV a ,b(x, y) - hxy, y) 

= (bx 2 -hx + a + b 2 ) 2 - 4(ax + b) (ax 2 + a 2 x + bx + ab). 

Then we have to compute the resultant of the above polynomials with respect to x. It 
always decomposes as the product of two quartic polynomials in h. Its expression is very 
large, so we only give it when a = 3 and 6=1/2. It writes as 
625 

(4/i 4 - 1176/i 3 + 308/i 2 + 287380/j + 1816975) pJh). 

65536 

It has four real roots, two for each polynomial. Some further work proves that the one that 
interests us is the smallest one of p^. g 
We also give an example when k = 3, but skipping all the details. 

Example 8. Consider the recurrence (pQ) with k = 3 and {a n } n as in (J3j) taking the values 
a = 1/2, 6 = 2 and c = 3. Then for any positive initial conditions x\ and Xi, V c \, a (%i, ^2) = 
h > V Ct b a (p) = h c — 15.9283. Moreover if the rotation number of F c ^, a associated to the 
oval {V C} b, a (x , y) = h} is irrational then the adherence of {x n } n is given by: 

• Three intervals when h £ (h c ,h*), where h* ~ 15.9614; 

• Two intervals when h £ [h*,h**), where h** ~ 16.0015; 

• One interval when h S [h**, 00). 

The values h* and h** are roots of two polynomials of degree 8 with integer coefficients that 
can be explicitly given. 
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5 Some properties of the rotation number function 



From Theorem [2] it is natural to introduce the rotation number function for F^ a : 

Pb,a ■ [h c ,oo) — ► (0,1), 

where h c := V& )(I (p), as the map that associates to each invariant oval {Vb >a (x,y) = h}, the 
rotation number Pb, a {h) of the function Fb a restricted to it. The following properties hold: 

(i) The function pb, a (h) is analytic for a > 0, b > 0, h > h c and it is continuous at h = h c . 
This can be proved from the tools introduced in Sec. 4]. 

(ii) The value pb, a (h c ) is given by the argument over 2ir of the eigenvalues (which have 
modulus one due to the integrability of i*&, a ) of the differential of Fb >a at p. 

(hi) p b ,a(h) = p a ,b{h). 

(iv) p a ,a{h) = 2p a (h) mod 1, where p a is the rotation numbej*] function associated to the 
classical Lyness map. Then, from the results of [3] we know that pn(h) = 3/5, that 
for a / 1, positive, p a ,a(h) is monotonous and lim/ ! ,_ i . 0O p a , a (h) = 3/5. 

Note that item (iii) follows because F a b is conjugated with Fb a via tp = F a which 
is a diffeomorphism of Q + , because tp~ 1 F aj bip = F~ F a FbF a = FbF a = Fb t a- Since V 
preserves the orientation, the rotation number functions of F a fi and Fb^ a restricted to the 
corresponding ovals must coincide. 

Similar results to the ones given above hold for F C ^ A an d its corresponding rotation 
number function. Some obvious differences are: 

Pc,b,a( h ) = Pb,a,c( h ) = Pa,c,b(h) Pa,a,a( h ) = 3 P( h ) mod 1) 

Pl,%A h ) = 2 / 5 i lim Pa,a,a( h ) = 2 / 5 - 

h^oo 

We are convinced that when a > and b > 0, 

lim Pb, a {h) = 3/5 and lim p c ,b,a(h) = 2/5, 
h— >oo h— >oo 

but we have not been able to prove these equalities. If they were true, by combining them 

with the values of the rotation number function at h = h c this would give a very useful 

information to decide if, apart from the trivial cases a = b = l(c = 1), there are other cases 

for which the rotation number function is constant. Notice that in these situations the 

* Notice that given a map of the circle there is an ambiguity between p and 1 — p when one considers its 
rotation number. So, while for us the rotation number of the classical Lyness map for a = 1 is 4/5, in other 
papers it is computed as 1/5. 
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maps Fb, a or F c h a would be globally periodic in Q~^~ • This information, together with the 
values at h c , also would be useful to know the regions where the corresponding functions 
could be increasing or decreasing. Finally notice that this rotation number at infinity is not 
continuous when we approach to a = or b = 0, where the recurrence and the first integral 
are also well defined on Q + . For instance po,o(p) = 2/3 and the numerical experiments of 
next subsection seem to indicate that for a > or b > 0, 



lim po >a (h) = lim p b ,o(h) 

h— >oo h— >oo 



5/8. 



Before proving Theorem [3] we introduce with an example the algorithm that we will 
use along this section to compute lower and upper bounds for the rotation number. We 
have implemented it in an algebraic manipulator. Notice also that when we apply it taking 
rational values of a and b and rational initial conditions, it can be used as a method to 
achieve proofs, see next example or the proof of Theorem [3j 

Fix a = 3, b = 2 and (xo,yo) = (1, 1). Then h = ^,3(1, 1) = 34. Compute for instance 
the 27 points of the orbit starting at (1, 1), 

. , . . . /9 17\ /89 788 

(xi,yi) = (4,6), (x 2 ,j/2) = I ) ' ( X3 ' m > = [ 54' 153 

and consider them as points on the oval {V2${x,y) = 34}, see Figure [U 




Figure 1: Oval of {V2^{x,y) = 34} with 27 iterates of i^,3- The label indicates the initial 
condition (1, 1), and the label k, k = 1, . . . , 26, corresponds with the fe-th point of the orbit. 

We already know that the restriction of i*2,3 to the given oval is conjugated to a rotation, 
with rotation number p := P2,3(34) that we want to estimate. This can be done by counting 
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0,582 




Figure 2: Lower and upper bounds for P2,3(34) obtained after following some points of the 
orbit starting at (1, 1). 



the number of turns that give the points (xj,yj), after fixing some orientation in the closed 
curve. We orientate the curve with the counterclockwise sense. So, for instance we know 
that the second point (x2>2/2) has given more that one turn and less than two, giving that 
1 < 2p < 2, and hence that p G (1/2,1). Doing the same reasoning with all the points 
computed we obtain, 



4 < 7p < 5 
8 < Up < 9 
10 < 19 p < 11 
14 < 26p < 15 



P G 

p e 
p e 
p g 



4 5 
7' 7 
8 9 
14' 14 
10 11 
19' 19 
14 15 
26' 26 



where we have only written the more relevant informacions obtained, which are given by 

the points of the orbit closer to the initial condition. So, we have shown that 

4 15 
0.5714 ~ - < ,02 3(34) < — ~ 0.5769. 
7 ' 26 

In Figure [5] we represent several successive lower and upper approximations obtained 

while the orbit is turning around the oval. We plot around six hundred steps, after skipping 

the first fifty ones. By taking 1000 points we get 

338 473 
0.5761246 ~ — < ,92,3(34) < — ~ 0.5761267, 



and after 3000 points, 

0.57612457 



338 , N 1472 

^<«,3(34)< — =.0.57612524. 



15 



In fact when we say that /»2,3(34) € (piow 5 Pupp)> the value p\ ow is the upper lower bound 
obtained by following all the considered points of the orbit, and p upp is the lowest upper 
bound. Notice that taking 1000 or 3000 points we have obtained the same lower bound for 
P2,s(34). 

Let us prove Theorem [3] by using the above approach. 

Proof of Theorem^ Consider a = 1/2, b = 3/2 and the three points 

i _/149 173 \ 2_(3 173 \ 3_( 1 173 

p "vioo'iooy' p vlo'ioo/ P ~ V iooo 'Too 

Notice that 

, n 10655559 
h\ := Kj/2i/ 2 (p ) = ^ 8.27, 

i/2,X/2\f ) 1288850 

, o, 9328327 
h 2 := W(P 2 ) = ^5- -44.93, 

, 1056238343 

h 3 := V 3/2 i/2 P = : - 3052.71. 

J/Ai/^H ) 346OOO 

Hence h c < h\ < h% < h%. By applying the algorithm described above, using 100 points of 
each orbit starting at each p*,j = 1,2,3, we obtain that 

/ , \ , /3 59\ , , /56 53 

P3/2,l/2(^l),P3/2,l/2("3j G ( 5'^ ) and ^3/2,1/2^2) € I — , — 

Since 59/98 < 56/93 we have proved that the function p 3 /2,i/2(h) has at least a local 
maximum in (hi,hs). Prom the continuity of the rotation number function, with respect 
a, b and h, we notice that this result also holds for all values of a and b in a neighborhood 
of a = 1/2, b = 3/2. | 

We believe that with the same method it can be proved that a similar result to the one 
given in Theorem [3] holds for some maps F c ^ a , but we have decided do not perform this 
study. 

5.1 Some numerical explorations for k = 2. 

We start by studying with more detail the rotation number function P3/2, 1/2(^)1 that we have 
considered to prove Theorem^ In this case the fixed point is p ~ (1.493363282, 1.730133891) 
and he = V bja (p) = 8.267483381. Moreover p b>a (h c ) ~ 0.6006847931. By applying our al- 
gorithm for approximating the rotation number, with 5000 points on each orbit, we obtain 
the results presented in Table 1. In Figure [3] we also plot the upper and lower bounds of 
Ps/2,i/2(h) that we have obtained by using a wide range of values of h. 
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Init. cond. (x, y) 


Energy level h 


Pkmih) 


PuppO) 


X 


h c ~ 8.2675 




~ 0.6006848 




~ 0.6006848 


1.3 


8.3068 


173 
288 


~ 0.6006944 


2938 
4891 


~ 0.6006951 


0.75 


9.2747 


1435 
2388 


~ 0.6009213 


2087 
3473 


~ 0.6009214 


0.3 


14.7566 


1548 
2573 


~ 0.6016323 


2285 
3798 


~ 0.6016324 


0.075 


44.9347 


657 
1091 


~ 0.6021998 


2354 
3909 


~ 0.6022001 


0.001 


3052.75 


2927 
4867 


~ 0.6013972 


86 
143 


~ 0.6013986 


5 • 10~ 6 


609716.07 


1832 
3049 


~ 0.6008527 


1409 
2345 


~ 0.6008529 


5 • lO" 256 


6.097 • 10 255 




| = 0.6 


2999 
4998 


~ 0.6000400 



Table 1: Lower and upper bounds of the rotation number P3/2,i/2{h), for some orbits of 
^3/2,1/2 starting at where p = (x,y). 




Figure 3: Lower and upper bounds for P3/2,i/2(M- O n the horizontal axis we represent 
— log 10 (/i) and, on the vertical axis, the value of the rotation number. Notice that for 
values of — log 10 (/i) smaller that 70 both values are indistinguishable in the Figure. 

For other values of a and b we obtain different behaviors. All the experiments are 
performed by starting at the fixed point p = (x,y), and increasing the energy level by 
taking initial conditions of the form (x,y), by decreasing x to 0. With this process we take 
orbits approaching to the boundary of Q + , that is lying on level sets of V^^ with increasing 
energy. The step in the decrease of x (and therefore in the increase of h) is not uniform, 
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and it has been manually tuned making it smaller in those regions where a possible non 
monotonous behavior could appear. 

Consider the set of parameters F = {(a, b), G [0, oo) 2 }, where notice that we also consider 
the boundaries a = or b = 0, where the map i^,o is we U defined. We already know that 
the rotation number function behaves equal at (a, b) and (b, a). Moreover we know perfectly 
its behavior on the diagonal (a, a) (when a < 1 it is monotonous decreasing and when 
a > 1 it is monotonous increasing) and that pi,i(h) = 4/5 and po,o(h) = 2/3. Hence a good 
strategy for an numerical exploration can be to produce sequences of experiments using our 
algorithm by fixing some a > and varying b. For instance we obtain: 

• Case a = 1/2. For all the values of b > considered, the rotation number function 
seems to tend to 3/5 when h goes to infinity. Moreover it seems 

- monotone decreasing for b G {1/4, 1}; 

— to have a unique maximum when b G {7/5,3/2}; 

— monotone increasing for b G {2,3}. 

• Case a = 0. For all the values of b > considered, the rotation number function seems 
to tend to 5/8 when h goes to infinity. Moreover it seems 

— monotone decreasing for b G {1/10, 3/10, 1/2}; 

- to have a unique maximum when &G {7/10, 3/4}; 

- monotone decreasing for b G {1, 5}. 

The above results, together with some other experiments for other values of a and b, not 
detailed in this paper, indicate the existence of a subset of positive measure in T where the 
corresponding rotation number functions seem to present an unique maximum. This subset 
probably separates two other subsets of T, one where Pb,a 

(h) is monotonically decreasing 
to 3/5 , and another one where pb, a {h) increases monotonically to the same value. The 
"oscillatory subset" seems to shrink to (a, b) = (1, 1) when it approaches to the line a = b 
and seems to finish in one interval on each of the borders {a = 0} and {b = 0}. Further 
analysis must be done in this direction in order to have a more accurate knowledge of the 
bifurcation diagram associated to the behavior of pb >a on T. 
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